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Abstract 

We introduce an integrable model for two coupled BCS systems 
through a solution of the Yang-Baxter equation associated with the Lie 
algebra sw(4). By employing the algebraic Bethe ansatz, we determine 
the exact solution for the energy spectrum. An asymptotic analysis is 
conducted to determine the leading terms in the ground state energy, 
the gap and some one point correlation functions at zero temperature. 
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1 Introduction 



The role of the Yang-Baxter equation in the study of quantum mechanical 
models has a long and distinguished history. Notable examples are the XYZ 
chain |]J , the t — J model at supersymmetric coupling [0 and the Hubbard 
model ||, each of which is both integrable and exactly solvable as a result 
of the formulation for each model through the Quantum Inverse Scattering 
Method (QISM). The key concept of the QISM is the notion of mutually com- 
muting transfer matrices, the existence of which is a result of the Yang-Baxter 
equation. In each of the above examples the Hamiltonian of the model is de- 
fined as the logarithmic derivative of the transfer matrix, and by the nature 
of the construction this yields a model defined on a one-dimensional lattice 
with nearest neighbour interactions, where it is possible for both integrability 
and solvability to hold for a variety of boundary conditions. 

The application of the QISM however can be applied on a much more gen- 
eral level. In the context of the present work, it is appropriate to mention, 
for example, the work of Gaudin [|J in relation to the construction of systems 
with long range interactions. Very closely related to Gaudin's Hamiltonians 
is the BCS model, the exact solution of which was quite remarkably found 
in 1963 by Richardson ||, while integrability was established much later by 
Cambiaggio et al. in 1997 |J. The exact solution of the BCS model has 
come under close scrutiny in recent years due to its application in the theory 
of metallic nanograins J7|. Specifically, the experiments of Ralph, Black and 
Tinkham || have shown that it is not valid to apply the BCS mean field 
theory for systems of nanoscale size. As a result, one has to turn to the 
exact solution of Richardson in order to conduct a reliable analysis. How- 
ever, the approaches adopted in || || make no reference to the Yang-Baxter 
equation or the QISM, (indeed QISM was not developed until many years 
after Richardson's work), and since historically the two facets of integrability 
and solvability have been intimately linked in the QISM framework, it was 
a natural question to ask whether the BCS model could be recast through 



this technique. An affirmative answer was given in |5], [K| with the surpris- 
ing result that the Z?-matrix solution of the Yang-Baxter equation which is 
needed in the construction of the BCS Hamiltonian is one of the simplest 
known examples, being that associated with the su(2) algebra. Given that 
a great volume of literature exists devoted to solutions of the Yang-Baxter 
equation associated with representations of simple Lie algebras, there is a 
vast opportunity to investigate generalised models. An important step to- 
wards this has already been achieved in [111] where a connection has been 
established between Chern-Simons theory and integrability of models associ- 
ated with an arbitrary Lie algebra, which is achieved through the Knizhnik- 
Zamolodchikov-Bernard equations. 

Such generalised models can be interpreted as coupled BCS systems, at 
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least in the sense that every simple Lie algebra can be generated by a system 
of simple roots which each form an su(2) subalgebra. An example of this was 
given in |L2| where the Lie algebra employed was so(5). In this instance, the 



model constructed describes proton-proton and neutron-neutron pairing as 
well as a coupling term for the scattering of proton-neutron pairs. Here we 
shall introduce a model based on the sw(4) Lie algebra symmetry which can 
also be interpreted as a nuclear system where there are now different types 
of pairing interactions. The Hamiltonian takes the form of two BCS systems 
which individually describe pairing interactions for the protons and neutrons 
and the scattering of bound proton pairs-neutron pairs, which is in contrast 



to the proton- neutron pairs of |12j . This interpretation is possible because the 
number operator for each system provides a good quantum number; i.e. the 
number operators are conserved. Therefore we can identify each BCS system 
with a particular distinguishable particle, which in this case are the protons 
and neutrons. It is worth remarking that this situation is inherently different 



to the pairing models described in [14] based on higher spin representations 



of the su(2) algebra. In these cases, the only good quantum number is 
the total number of particles in the combined system. There, individual 
particle numbers are not conserved and thus the models can be interpreted 
as describing a Josephson tunneling phenomena. 

The paper is organized as follows. In Section 2 we present the construction 
of the model through the QISM. In section 3 the exact solution of the model is 
given by means of the algebraic Bethe ansatz. An analysis of the asymptotic 
solutions of the Bethe ansatz equations is presented in Section 4 , where the 
ground state energy, the gap in the spectrum of elementary excitations, as 
well as the derivation of some correlation functions in this asymptotic regime, 
are presented. A summary of the main results can be found in Section 5. 

2 Coupled pairing Hamiltonian and integra- 
bility 

Let us begin by introducing the following Hamiltonian 

n 



H = BCS(l)+BCS(2)-gJ2K( 1 ) b t( 2 M 2 M 1 ) 
n 

+ 5 ^&+(l)& fc (l)( % (2)-n fe (2)) 2 



n 

+gJ2b+(2)b k (2)(n J (l)-n k (2)f, (1) 
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where 

n n 
£CS(a) = J2 2e j n j( a ) ~9 h t( a ) h k{a)- (2) 

Above the operators bj(a), b~j"(a) are the annihilation and creation operators 
for the hard-core bosons (or Cooper pairs) in system a, and j refers to the 
single particle energy level with energy e^. We will assume that the values 6j 
are distinct. Further, g is a coupling strength constant for the scattering of 
Cooper pairs and rij(a) = b^(a)bj(a), is the Cooper pair number operator. 
As in the case of the usual BCS system there is a blocking effect (e.g. see 
0), as there is no scattering of any unpaired states. For each level j there 
are actually sixteen local states, but the nature of the Hamiltonian means 
that only on a subspace spanned by four of these states, where there are no 
unpaired states, is the scattering non-trivial (see (|7|) below). Hereafter we 
will restrict our analysis to this subspace. 

On this restricted subspace the operators b1-{a) = c^(a)c^(a), bj(a) = 

Cji(a)cj^(a), where c Jcr , c] CT , & =t, I, are the familiar fermion operators, sat- 
isfy the hard-core boson relations 

(6+(a)) 2 = 0, [bj(a), b+(b)] = 5 ab 6 jk (l - 2bj(a)b j (a)), 
[6,(a),6 fc (6)] = [6+(a),6+(fe)]=0, for k^j. 

We can see from the Hamiltonian expression that the exchange interaction 
of Cooper pairs in one system depends on the number of Cooper pairs in the 
other system. For example, if in system (1), the level j is empty and the level 
k is occupied by one Cooper pair, just for certain configurations of system 
(2) it is possible that this Cooper pair in (1) scatters from level k to j. This 
means that the Hamiltonian ([!]) presents naturally some "selection rules" for 
the scattering of states. We illustrate these configurations to indicate the 
possible pair scatterings 
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In addition, the double-pair scattering terms of the form 



(2) 




are also present. What the above indicates is that besides the number of 
Cooper pairs being conserved in each system, the number of double pairs 
(to be more precise, the number of energy levels which are completely filled) 
is also conserved. This can be seen in each of the scattering processes de- 
picted graphically above. In each case the scattering does not overall change 
the number of completely filled levels. There are further symmetries in the 
Hamiltonian. For example, there is a reflection symmetry which interchanges 
the labels 1 and 2 for the two BCS systems. This arises as a result of a global 
so(3) © symmetry that the model possesses, which will be made more 
clear later. In that which follows we shall first discuss the integrability of the 
Hamiltonian ([[]) in the context of the QISM. 

In order to built up a mechanism to construct an integrable su(A) pair- 
ing model, let us first recall the quantum i?-matrix associated with the Lie 
algebra su(4), which acts in the tensor product of two 4-dimensional spaces 
V ® V and can be written as 



R(X) 



(\.I®I + r)P) 
(A + v) 



(3) 



Above A is the usual spectral parameter, P is the permutation operator 



with matrix elements P. 



<5 q ,5<Vp a )/^;7;^ = 1)2,3,4 and rj is the 



quasiclassical limit parameter; i.e. 

lim R(X) = I <g> I. 

?7— »o 

It is known that this i?-matrix satisfies the Yang-Baxter equation (YBE) 

R 12 (X - fi)R 13 (X)R 2 3^) = R 2 3(v)Ri3(X)Ri2 (A - fi). (4) 

The R-matrix may be viewed as the structural constants for the Yang-Baxter 
algebra generated by the monodromy matrix T(A), namely, 



Ru(X -fi)T (A) T (fi) =T {fi) T (X)R 12 (X - fi). 



(5) 
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Consequently, the f?-matrix (^) allows us to construct a realization of the 
monodromy matrix through 

T(A) = G i?oc(A - e n ) ■ ■ ■ Goi? i(A - e x ). (6) 

Here the subscript denotes the auxiliary space and G satisfying 

[R, G <g> G] = 

is a class of c- valued solutions of the YBE (^) . Asa consequence of the Yang- 
Baxter algebra (0), the transfer matrices t(X) = tr T(A) mutually commute 
for different values of the spectral parameter A. This transfer matrix is the 
starting point in the construction of a su(4)-type Gaudin Hamiltonian, from 
which we can obtain the su(4) pairing Hamiltonian, as will be shown below. 
For this purpose we make the following identification for the basis states 



|1>= |0> 



|2> = b+(l)&+(2)|0>= 4f4^ 



|3>= & + (l)|0) 



4f 



|4)= 6+(2)|0> 
and choose the G-matrix to be given by 



G = exp 



2r}{l - n(l) - n(2)) 



/ exp'' — 



tog' 



4t 



\ 



exp(^) 





V o 



Qg 

1 



1/ 

to construct the transfer matrix t(X). It can be verified that 

' ' " ( 'j h'.i-j- -('.i ~ e i+i) G r 
Next, taking the quasiclassical limit, we find 

R hk {X)\ v ^ = I ®I + r 1 r hk {X) + 0{rf) 

G;Uo = / + ^(l-n J (l)-n,(2)) + 0(r / 2 ) 



(7) 



(8) 



(9) 

(10) 
(11) 
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where ^(A) 



-. Thus it follows that 



t(ej)\ v ^ = 1 + 7? 



n 

£ t4 



fc = i 



+ 



(12) 



where 



r, = -(l-n J (l)-n,(2))+ £ ^ g ' 



4 p/3a 
fc 



k = 1 



e.7 - e fe 



(13) 



Here i? Q/3 = |a)(/3|, a,/3 = l,---,4 are the Hubbard operators. An imme- 
diate consequence from the Yang-Baxter algebra (|5|) is that [Tj,r fc ] = 0. In 
addition, as a result of the so(3) © u(l) symmetry mentioned earlier, it can 
be shown that there are extra conserved operators K and % such that 

[rj, K] = [ Tj , X ] = [K, X ] = 0. 

Above, K is the Casimir operator of an so(3) subalgebra acting on the f2-fold 
tensor product 



j,k v J 



(14) 



where (L°, L + , L ) are the basis elements of this canonical so(3) subalgebra 



L+ = E u = b + (l)b(2), 
IT = E 43 = 6+(2)&(l), 
L° = E 33 -E u = n{l)-n{2). 



(15) 



The -u(l) operator x explicitly reads 

X>f + 

n 

^(n,(l)-n,(2)f 



X 



(16) 



Any Hamiltonian which is defined in terms of the mutually commuting 
set of operators 

{Tj, K, X } (17) 
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will necessarily be integrable where the operators in ([T7p represent the con- 
stants of the motion. By making the following choice 

£7 o £7 o £7 

+| X (x - ft) + 2 ^ £i + ^- - 2^ (18) 

3 

we produce the Hamiltonian (ffj). In order to determine the energy spectrum 
of this model, we will need to determine the eigenvalues of the conserved 
operators, to which we turn next. 



3 Bethe ansatz solutions 

Besides proving the integrability of the model, we can also obtain its exact 
solution from the algebraic Bethe ansatz for the standard sit (4) vertex model 
constructed from the R- matrix (|^). Employing the nested algebraic Bethe 
ansatz jn| we can obtain the eigenvalue of the transfer matrix (|6]) as 

N 

AMA,}K}) = e» II * 

i=i u Ui 

Q, n N M n 

+e -f f[ v Z e A ~ I TT v ~ V A ± 1 TT v ~ Ul ~ 2 



n 



91 — 4- — ■*- ■*- it — in j 4- n 



. V - €i + j ±A ; V-Wj+rj 
i=l * j=l J 

Above the parameters Vj , 77 m and Wk satisfy the Bethe ansatz equations 



n „ N M 



n 



e -^JTV j -e i - 1 jrV j -V l + V jrV j -U l - 1 

J. J. 71. p. _|_ _i J. J. «). 1), — n J- J- >?!■ — >?/, _l_ _i 



± TJj - Cj + \ A-J TJj - V« - 77 A. 1 ^ - M/ 

at . v M Q 



.£2 T-T «m - «i + | _ tt «m - i*i + »? TT U m ~ W i " 2 

f-J-^m-TJi-f 11 Mm - 7j- - 77 11 Wm - lyj + 2 

1=1 z 1=1 t=l z 

M ri Q 

y-r Wfc - Ui + j _ y-r W k -Wi + 1] 

f=l W k-Ul~l ^W k - Wi - T] 

Ijtk 
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j = l,---,N, m = l,...,M, fc = l,...,Q. 

Defining iV(a) = 2\^ =1 n.,(a), we can readily determine that the quantum 
numbers N, M and Q are given by 

N = iV(l)+JV(2)-JV(l)JV(2), 
M = JV(l)+JV(2)-2JV(l)iV(2), 
Q = iV(2)-iV(l)JV(2). (20) 

The eigenvalues of the integrals of motion Tj (|T3| ) can be obtained from 
the expansion of the eigenvalue of the transfer matrix (|l^) in the parameter 
r). Explicitly, the eigenvalues of Tj are given by 



9 N n 



at n 

E^t+Et^ < 21 > 



O c — ' Vi — €j * — ' 6a — 6u 

^ 1=1 1 J k=l 3 

where the parameters satisfy the following equations 

n m n 



~+y — + y — =2y — . 

9 ^Vj-6i ^Vj-Ui ^Vj-Vi 



1=1 



N M 



9 Urn-Vi *-i Ui-Um f^f U m ~ lUj 

1=1 1=1 1=1 

M Q . 

E^ = 2 E^— > ( 22 ) 

j = l,---,N, m = l,---,M, fc = l,...,g 



We will also need the eigenvalues of the operators K and x- Through use of 

( has eigenvalue M while 

-{M -2Q)(M -2Q + 2). 



lq , |20|) we find that x has eigenvalue M while the eigenvalues of K are 
1 
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Finally, utilizing ( 21"D and noting the following identities which can be 
derived from (|22|) : 

M 



M N 

EE 



k =i i=i w «~ u i t^il=1 Wk - Wl 
" w A ' 2M 



m=l 4=1 



Q M Q M 

M N M N 

m=l 1=1 m=l t=l 

N 0, NO. 

^ 6i - Va ^— ' ^ £i - Va 

j=l i=l J j=l i=l J 

Q M 

k=l 1=1 K 1 

M M N M Q 

m=l m=l i=l m=l 1=1 

. N N Q N M 

we can present from the relation ([Tj|) the eigenvalue of the Hamiltonian ([!]) 

as 

JV M 

E = ^Y, Vi -2^« m -#(2N-3M). (23) 

i=l m=l 

Let us make some small remarks about the degeneracies of the spectrum. 
Though the eigenstates of the Hamiltonian have not been made explicit here, 
it can be deduced by the standard arguments (e.g. p6| ) that each is a highest 



weight state with respect to the so(3) symmetry algebra (|15D . In particular, 
the highest weight which is given by the eigenvalue of the operator L° is 
M — 2Q, so we can conclude that the multiplet generated by ( |15D acting 
on this highest weight state has dimension M — 2Q + 1. Therefore for each 



solution of (22) with given N, M and Q, the corresponding energy level has 



degeneracy M — 2Q + 1. 



4 Asymptotic solutions 



As the Bethe ansatz equations (p2[) take the form of coupled non-linear equa- 
tions it is unlikely to find analytic solutions, and one tends to resort to nu- 
merical analysis. It is however possible to conduct an asymptotic analysis 
for small values of the coupling parameter g. Below we undertake this for 
the ground state of the system and some elementary excitations. 
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For the ground state, we consider first the case g = 0. Letting N% and 
N 2 denote the number of Cooper pairs in each system, then it is clear that 
the ground state corresponds to filling the Fermi sea, which is illustrated in 
Fig. (A) below, where without loss of generality we assume that N\ > N 2 . 
For small g 7^ we see that the ground state will be described by a solution 
of (P|) with N = N u M = JVi - N 2 , Q = 0. 




Fig. (A) 



Thus the Bethe equations ( p2|) reduce to two levels for the parameters Vj and 
u m . For a small g > it is appropriate to consider the asymptotic solution 



Vj = ej+gdj + g 2 ^, j = l,---,N 1 , (24) 

u m = ejVi-m+i + 9 a m + g 2 Pm, m = 1, ■ ■ ■ , N x - N 2 . (25) 

Substituting these into the Bethe equations (|22| ) with the configuration N = 
N u M = Aq - N 2 , Q = 0, one can find that 



~ 9,9 
Vj pa e,- 

3 3 4 16 



9 , 9 

v 3 ~ e J'~2 + T 



n ^ n 2 ^ 



£i — e< — €j — £j 



i= AT 1 + l 3 1 J =1 



, 3 < N 2 , 



a Nl 1 



(=^2 + 1 



e 7 - e / 



A'l 



A'2 



E — - — + E 



(26) 



, 3 > N 2 , (27) 



1 



e iVi-m+l — 



- £ 

Z=iV 2 +l 



m — 1, ■ • ■ , Ni — N 2 



ii 



The asymptotic ground state energy is deduced to be given by 



E « 4^ £j -2 e l -(N 1 + 2N 2 )g 

j=l l=N 2 +l 



+ 



9_ 

4 



n 2 n 



/ ^— ✓ e . _ e . ^~ ^— ✓ 1 ^— ✓ e . _ e . 

j=l i=JVi+l J * J=iVa+l \i=N 1 +l 3 % 



N 2 



E" 



• 1 e i _ e * 

8=1 J 



It is important to point out that from the above ground state energy we can 
infer some results about the asymptotic behaviour of zero temperature corre- 
lation functions. Specifically, by employing the Hellmann-Feynman theorem 
we have that 



1 f)W 

Wl) +ni (2)> = -^. 



The result obtained is 



(n i (l) + n i (2)) 
(ni(l)+7ii(2)) 
(n i (l)+n i (2)) 



.2 / ^2 



E 



,e,- — e 

j=JVi+l v J 



Vi 

E 

j=V 2 +l 

v 2 

-E 



e,- - Es 



e-i - e, : 



2 - 



T~~\2 + XI 



Vj=JVi+l ^' ' 



j= N 2 +l & €i 



for % > Ni, 
for iVi > % > N 2 , 
for i < N 2 . 



Next let us consider a possible excitation which can be obtained by break- 
ing one Cooper pair in BCS(2). In the g = case, the excited state is depicted 
in Fig. (B). 

(1) (2) 




Fig.(B) 

For non-zero g, we choose N = Nx — 2, M = N\ — N 2 — 1 and block the levels 
with energy €n 2 , en 2 +i- From the asymptotic solutions (pop, (p7|) and (|28|). 
we obtain the excitation energy 



Ni 



E 

j=l l=N 2 +2 



ei - en 2 - en 2 +i ~ {Nt + 2N 2 - A)g 



12 



9 2 



■jv 2 -i n 

E E 

_ j=i i=JVi+l 



E E 

j=N 2 +2 \i=N 1 +l 



N 2 -l 

E 



, e i ~ €i 
■1=1 J 



Therefore the gap obtained through the breaking a Cooper pair in BCS(2) 
is given by 
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E 



U — CJV; 

2 



+ E - 



i=JVi+l 
iV 2 



+E 



e N 2 +l 



Another possibility for an excitation at g = is to break the Cooper pair 
at level in BCS(l) (see Fig.(C)). For g 7^ 0, the configuration should be 
accommodated as N = Ni — 1, M = Ni — N 2 — 1 with the levels ejv 15 e7v 1+ i 
blocked. 




Fig.(C) 

We find the gap obtained in this case to be given by 



A, 



9 

— € N! + 9 + — 



' N 2 

E 



Ni-1 

+ E 

j=JV 2 +l 



ejvi+i 



N 2 

+ E 



£TVi+1 



+ Er^ 



j=JVi+l 



For the asymptotics to be valid we require from (|(| |27|, p8|) that g << 
for all i. However in the above two expressions for the energy gaps the value 
of g may still be larger than the spacings e^+i — 6i of the single particle 
energies. So we can see from the analysis when this is the case that the 
inclusion of the pairing interactions does produce a significant gap between 
the ground state and the first excited states. 
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5 Conclusion 



To summarize, we have constructed an integrable pairing Hamiltonian based 
on the su(4) Lie algebra. This model can be interpreted as describing two 
coupled BCS systems of different types, such as for protons and neutrons 
in a nuclear system. The Bethe ansatz equations and the energies of the 
model have been calculated. For small values of the coupling parameter g, 
we asymptotically analyzed the ground state and elementary excitations, and 
the expectation values for the occupation numbers. An open problem that 
we will address in the future is the exact calculation of form factors and 
correlation functions using the techniques developed by Babujian et al. in 
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